Bubbling in delay-coupled lasers 
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Figure 1: Schematic setup. 
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We study chaos synchronization of two lasers which are delay-coupled via an active or a passive 
relay. While the lasers are synchronized, their dynamics is identical to a single laser with delayed 
feedback for a passive relay and identical to two delay-coupled lasers for an active relay. Depending 
on the coupling parameters the system exhibits bubbling, i.e., noise-induced desynchronization, or 
on-off intermittency. We associate the desynchronization dynamics in the coherence collapse and 
low frequency fluctuation regimes with the transverse instability of some of the compound cavity's 
antimodes. Finally, we demonstrate how, by using an active relay, bubbling can be suppressed. 



Synchronization phenomena of coupled nonlinear os- 
cillators are omnipresent and play an important role in 
physical, chemical and biological systems [H, @|. Under- 
standing the synchronization mechanisms is crucial for 
many practical applications. One of the most interest- 
ing and challenging phenomena when coupling nonlinear 
systems is the synchronization of chaotic dynamics 
In order to characterize the synchronization effects, sta- 
bility properties are a key issue. Noise can, for instance, 
cause intermittent desynchronization. This behavior is 
called bubbling [3] and has been observed for example in 
optical [US] and electrical [3] systems. 

Semiconductor lasers are of particular interest in 
the study of chaos synchronization. The synchroniza- 
tion properties may facilitate new secure communication 
schemes. However, if two identical semiconductor lasers 
are optically coupled over a finite distance, it has been 
observed that the coupling delay leads to spontaneous 
symmetry breaking, and only generalized synchroniza- 
tion of leader-laggard type occurs 0]. A passive relay 
in form of a semitransparent mirror or an active relay 
in form of a third laser in between the two lasers have 
been shown to stabilize the isochronous synchronization 
solution rendering such configurations at- 

tractive for chaos encrypted communication and other 
chaos based applications fl3| . 

In this Letter we show that bubbling and on-off in- 
termittency occur in both relay setups. In the coher- 
ence collapse (CC) and in the low frequency fluctua- 
tion (LFF) regime, we find that bubbling is caused by 
transversally unstable external cavity modes (ECMs). In 
the LFF regime the localization of the transversally un- 
stable modes in the synchronization manifold (SM) re- 
sults in desynchronization during power dropouts, which 
has also been observed in unidirectionally coupled lasers 
[3]. For the active relay we find that bubbling can be 
suppressed by stronger pumping of the relay laser. 

We consider two identical systems which are delay- 
coupled via a relay (Fig. [I]). The relay may be an active 
element or a passive element which merely distributes 
the arriving signals between the systems. Each system 
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G R n are the state vectors of the system 
j and the relay, respectively, f is a nonlinear function, 
K is the relay-to-system coupling matrix, and r is the 
propagation delay between system 1 and system 2. The 
overdot denotes the derivative with respect to time t. 
For the active relay we consider the equation 

Y = g(Y) + |LXi(t - r/2) + ±LX 2 (i - r/2), (2) 

where L is the system-to-relay coupling matrix and the 
function g describes the internal dynamics of the relay. 
For the passive relay we consider the algebraic equation 



Y(t) = ipti(i - r/2) + X 2 (t - r/2)]. 



(3) 



Equation {1} together with the relay equation §2$) or J3]) 
allow for an isochronous (or zero-lag) solution Xi(f) = 
X2(£), respectively. The SM is thus invariant. To analyse 
the stability of this solution we introduce a symmetric 
variable S = ^(Xi + X2) and an antisymmetric variable 
A = i(Xr 



v „i.i — X2). Equation ([I]) can then be rewritten in 
the new variables 



S = \ [f (S + A) + f(S - A)] + KY(t - r/2), (4) 
A= i [f (S + A) - f (S - A)] . (5) 

Note that due to the symmetric coupling the delay terms 
and all the coupling parameters in Eq. J5]) vanish. Equa- 
tion J5]) taken at A = has a solution A = which 
represents the isochronously synchronized state. Its sta- 
bility is determined by linearizing Eqs. (@]) and |(5]) in the 



variable A around A 
to the SM: 



0, i.e., we linearize orthogonal 



S=f(S) + ifY(t-r/2), 
A = Df (S)A. 



(C) 
(7) 



Here, Df(S) denotes the Jacobian of f evaluated at the 
position S. Since S depends on time, Eq. (7]) constitutes 
a time-dependent variational equation. 

For both relay types the dynamics within the SM re- 
sembles the dynamics of a single system with either self- 
feedback (passive relay) 



S = f (S) + KS(t - t) 



or coupling to the active relay 



S = f(S) + A^Y(<-r/2), 
Y = g(Y)+LS(t-r/2). 



(8) 



(9) 



In both cases the stability of the synchronized solution 
is governed by Eq. ((7|)- However, the trajectory S(t) will 
be different and the synchronized state may thus have 
different stability properties. 

Bubbling occurs 0,113] when an invariant set /, for ex- 
ample a periodic orbit, in the SM is transversally unsta- 
ble, while the chaotic attractor in the SM is still transver- 
sally stable, i. e. the largest transversal Lyapunov expo- 
nent of the attractor is negative, Aj_ < 0. In this situa- 
tion the trajectory can be pushed towards the unstable 
set by noise and leave the SM. If there is no other at- 
tractor present, the trajectory will eventually come back 
to the SM and the systems will synchronize again. The 
point where the invariant set I loses its transverse stabil- 
ity is called bubbling bifurcation, while the point where 
the attractor itself becomes unstable is called blow-out 
bifurcation. 

For semiconductor lasers the dynamics of each system 
is governed by the dimensionless Lang-Kobayashi rate 
equations 3, 13] 



Ej = |(1 + ia)n 3 E, + Ke lip E v (t - r/2) 



(10) 



p-n j -(l + n i )|.E J f 



Here, Ej and £V are the complex electric field ampli- 
tudes of the jth system and the relay, respectively, rij is 
the excess carrier density, a is the linewidth enhancement 
factor, p is the pump current, and the timescale param- 
eter T = t c /t p is the ratio of the carrier (r c ) and the 
photon (t p ) lifetime. For simplicity we choose the feed- 
back phase tp = 0. Note that in general one could also 
include coupling phases in Eq. (J3j) • This leads to interfer- 
ence conditions of all phases which have to be satisfied 
for isochronous synchronization. In our simulations we 
consider the spontaneous emission noise via a complex 
Gaussian white random variable Fe(£) with the covari- 
ance (F E (t) F E (t')*} = f3{n + n )S(t -t'), where n = 10 




is the carrier density at threshold and {3 = 10 -5 is the 
spontaneous emission factor. Carrier noise has not been 
taken into account at this level. 

If the relay is realized through a semitransparent mir- 
ror (passive relay), the dynamics within the SM is given 
by Eqs. QEJ) with Ey(t - r/2) = E 3 (t - r), i.e., an 
effectively decoupled laser. For this configuration we 
calculate the maximum parallel Lyapunov exponents A|| 
(within the SM) as well as the maximum transversal Lya- 
punov exponents A_l by simulating the dynamics in the 
SM without noise and applying the method developed 
in [3]. Figure [2^ displays the Lyapunov exponents as 
a function of the feedback strength K. There are two 
blow-out bifurcations at K w 0.008 (Bl) and at 
K w 0.09 (-B2), where Aj^ changes sign and the chaotic 
attractor loses its transversal stability. Similar behav- 
ior is found for an active relay (Fig. [2}}) . Over a wide 
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Figure 2: (Color online) Maximum transversal Lyapunov ex- 
ponent Ax (red dashed) and maximum parallel Lyapunov ex- 
ponent A| (blue solid) as a function of the feedback strength 
K for a) passive relay b) active relay (p re iay = 4.0). At the 
two blow-out bifurcations Bl and Bl the maximum transver- 
sal Lyapunov exponent of the chaotic attractor changes sign. 
Other parameters: T = 200, p = 1.0, r = 1000, a = 4, tp = 

range of K (Fig. Eh) in which the attractor is stable and 
the dynamics is chaotic, we observe bubbling induced by 
spontaneous emission noise. In these regimes, when the 
noise is switched off in the simulations, the two lasers 
stay perfectly synchronized. In the regime with A^ > 
we observe desynchronization bursts even without noise, 
i.e., the system exhibits on-off intermittency. Figure [3K 
depicts the bubbling behavior for values of K above B2 
where the laser operates in the CC regime. Figure [3b 
corresponds to a lower pump current, where the synchro- 
nized lasers operate in the LFF regime. In this regime 
bubbling only takes place during the power dropouts. In 
both cases, when the noise amplitude is decreased, the 
desynchronization peaks occur less frequently, the maxi- 
mum height, however, does not decrease. 

We now relate the desynchronization dynamics to the 
transverse stability of the ECMs in the SM. These modes 
organize the dynamics in the SM in the CC and the LFF 
regime. The ECMs are rotating wave solutions of the 
form E(t) = Acxp(io;t) and n(t) — n with constant val- 
ues A, oj and n. They are well studied solutions of 
the Lang-Kobayashi equations and are located on an el- 
lipse in the (w, n)-plane (see inset of Fig.HK). The modes 
on the top and bottom half of the ellipse are called modes 



and antimodes, respectively. 




t/10 4 

Figure 3: Carrier density of the symmetric variable 
ns = 2(^1 + 712 ) an< ^ intensity difference |7i — h \/(Ii + I2) 
(normalized by the mean intensity) representing the devia- 
tion from the synchronized state vs. time, a): Bubbling in 
the coherence collapse regime (p =1.0). b): Bubbling in 
the low frequency fluctuation regime during power dropouts 
(p = 0.1). Other parameters: T = 200, K = 0.12, r = 1000, 
a = 4. 




[<t>s(t)-Mt-r)]/r [<f> s (t)-Mt-r)]/T 

Figure 4: (color online) Projection of the dynamics of the 
symmetrized solution ns, Es = A s exp{i<f>s) (black trajec- 
tory) onto the (w,n)-plane for a) passive relay (p = 1.0), 
b) active relay (p = 1.0, p re ia y — 4.0). Transversally stable 
(blue circles) and transversally unstable (red squares) ECMs 
are also shown, a) The competition between chaotic itiner- 
ancy and antimodes leads to bubbling during global antimode 
dynamics. Yellow diamonds mark the onset of desynchroniza- 
tion. Solid and dashed parts of the trajectory correspond to 
synchronized and desynchronized periods, respectively. The 
inset in a) shows the ECM ellipse and bubbling dynamics in a 
larger range, b) The system evolves around the transversally 
stable compound laser modes and bubbling is suppressed. Pa- 
rameters as in Fig. [3] 



The transverse stability of an ECM is governed by the 
variational equation J7]) where S(t) is the ECM solution. 
To determine the stability, we transform the laser equa- 
tions into a rotating frame Eq = Ecxp(— icut). In 
these coordinates, an ECM E = Aexp(iiot + iip) is trans- 
formed into a family of fixed points Eq = Aexp(itp). 



Splitting the complex electric field Eoj 



iyj and 



using the vector Xj = (xj,yj,n,j) Eqs. lfT0|) can be writ- 
ten in the form of Eq. {J} and the above analysis applies. 
In particular, the Jacobian in the rotating frame is given 

by 



Dt = 



4(i 



l — Ld 2 n 

+ n)2x -±(\ + n)2y 



zOLX 



\ay 



-±(l + x 2 +y 2 ) 



The eigenvalues of the Jacobian then determine the 
ECM's transverse stability. Figure displays the po- 
sition of the ECMs in the (w, n)-plane and their stability 
for a choice of parameters. The black trajectory displays 
the projection of the symmetric variable ns. 

The bubbling behavior in the CC regime and the corre- 
lation of the desynchronization with the power dropouts 
in the LFF regime can be understood as follows. In 
the CC regime the dynamics comprises chaotic itiner- 
ancy among the modes and global antimode dynamics 
(see Fig. [4^). The modes involved in the chaotic 
itinerancy are transversally stable (blue circles) . The an- 
timodes on the other hand are transversally unstable (red 
squares). Thus, when the trajectory approaches the anti- 
mode, noise can lead to desynchronization and bubbling 



occurs. The yellow diamonds in Fig. 4a mark the on- 
set of desynchronization, showing that bubbling always 
occurs in the vicinity of the antimodes (independent of 
the power). Please note that due to the role of noise 
not every approach to an antimode results in a bubbling 
excursion. 

In the LFF regime [13] the dynamics is similar. The 
intensity buildup process in between power dropouts is 
characterized by chaotic switching between different at- 
tractor ruins (ghosts) of unstable ECMs with a drift to- 
wards the ECM with minimal n. All ECMs involved in 
the buildup process are transversally stable and we ob- 
serve no desynchronization. After a transient time, a 
power dropout takes place. During the dropout the tra- 
jectory collides with an antimode in a crisis. Again, the 
vicinity to transversally unstable antimodes - rather than 
the drop in power - leads to bubbling behavior. 

The transverse stability of the ECMs depends on the 
laser and coupling parameters as well as on the param- 
eters of the particular ECM. Note that modes and an- 
timodes are not necessarily transversally stable or un- 
stable, respectively. The modes on the lower right-hand 
side in Fig. 0K, for instance, are transversally unstable. 
With decreasing coupling strength K, more modes be- 
come transversally unstable until the whole chaotic at- 
tractor loses its transversal stability. This leads to the 
blowout bifurcation B2 in Fig. [H 

For very large feedback strength K > 0.7 no bubbling 
occurs. This can be expected, since A_l becomes more 



negative with increasing K (Fig. [2]). Note, however, that 
the parallel Lyapunov exponent Am approaches zero with 
increasing K and the chaoticity decreases, making this 
situation less interesting for chaos-based applications. 

If the elements are coupled via an active relay, the 
synchronized lasers behave like two delay-coupled lasers 
(see Eq. ©). If we choose f = g and K = L, we ob- 
tain a system of two identical mutually coupled semicon- 
ductor lasers, which has been studied before (il. li-l Eil]. 
Such a system has rotating wave solutions of the form 
E s (t) = A s exp(iwi), Ey{t) = Ay exp(itut + iip), n s (t) = 
ns, n-y(t) ~ ny, called compound laser modes (CLMs). 
Their spectrum is more complex than for the ECMs: 
besides the synchronized solutions (which correspond to 
ECMs), there exist antisymmetric modes, for which the 
relay and the synchronized solution are in anti-phase 
(ijj = 7r), as well as asymmetric modes where the relay 
has a different intensity than the outer lasers. 

The positions of the transversally unstable modes are 
close to those of the ECMs of a single laser in the (w, n) 
parameter space. Also the dynamics of three identical 
coupled lasers is similar to the behavior in the presence 
of a passive relay. Indeed, we find bubbling in both the 
LFF and CC regime. 

In the experiments reported in [12] all the coupling pa- 
rameters in the setup are chosen identical, i.e., L = 2K 
in Eqs. (TJ and ([2]). But also in this case we observe 
qualitatively similar laser dynamics, with a trajectory in 
parameter space coming close to the transversally unsta- 
ble CLMs. 

To suppress the bubbling while maintaining strong 
chaos, we apply a sufficiently higher pump current to 
the relay laser (p re iay = 4.0) than to the outer lasers 
(p = 1.0). For this configuration we have calculated 
A|| w 0.026, Aj_ ~ —0.032, confirming that the system 
is in the chaotic regime (cf. Fig. Eb)- The system still 
itinerates among the compound laser modes, but there 
is no global antimode dynamics. Moreover, in contrast 
to the behavior for the symmetric case p Ye ] ay = 1.0, the 
active relay now suppresses the bubbling and there is no 
desynchronization (see Fig. Sb). Inspecting Fig. we 
can conclude that the CLMs involved in the dynamics are 
indeed transversally stable. If the middle laser is pumped 
less strongly than the outer ones, the opposite effect is 
observed. 

In conclusion, we have demonstrated a mechanism for 
desynchronization by bubbling in a very general setting 
of two delay-coupled lasers with either passive or active 
relay. We have shown that in the CC and LFF regimes 
the occurrence of bubbling is related to the transverse 
instability of some of the compound cavity's antimodes, 
and that, by tuning of the active relay, it is possible to 
suppress the bubbling. This may have significant impli- 
cations for chaos synchronization of lasers and its appli- 
cations. 
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